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Abstract 



We define quantum field theory by taking the Lagrangian action to be given as a sequence 
of mathematically well-defined functionals written in terms of operator fields fulfilling given 
local commutation relations. The renormalized solution fields have a fully defined Fock space 
expansion and are multi- local; thus Haag's theorem does not apply, i.e., the interaction picture 
exists. Also, the formalism allows immediately the definition of a wave function and the 
description of many-body bound-state systems. 
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In the present paper we will show that it is possible to formulate relativistic quantum field theory 

(QFT) in a mathematically rigorous way in that every step in the chain of reasoning, beginning 

with the Lagrangian and terminating with the S-matrix can be rigorously defined and that all 
Q-i. 

the quantities participating in this chain exist as fully specified mathematical objects. In this 
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formulation the interacting theory joins smoothly to the non-interacting case. At the same time 
the treatment of bound composite system presents no difficulties; in particular, it allows the def- 



inition of a system wave function which is fully analogous to the wave function of non-relativistic 
quantum mechanics. In detail, the modifications from the conventional formulation of QFT are: 
(i) the basic fields have to be defined as local quantum fields; (ii) the Lagrangian functional 
must be replaced by a generalized functional (short: g- functional) , which is the analogue of the 
replacement of functions by generalized functions (short: g-functions), see below and Appendix A. 

Of the above mentioned two modifications, only the first is new; the second has been implied 
but only rarely used explicitly. Conventionally it is postulated that the "renormalized" fields 
are local fields. In the present formulation the basic fields of the Lagrangian, i.e., the so-called 
"free" fields, are defined to be local fields, and to obey once and for all fixed commutation 
relations. These fields are used to construct a Fock space, and all quantities of the theory are 
expanded in that basis. This way the "renormalized" fields, i.e., the solution fields, as the result 
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of the interactions acquire structure, i.e., they are non-local, and are not unitarily related to the 
"free" fields. Hence Haag's theorem |IJ does not apply in our formulation, i.e., the interaction 
picture exists. This way it is this first of the modifications which is the basis for the resolution 
of the contradictions inherent in the conventional formulation. One may say that the present 
re-formulation entails the reversal of the roles of the "bare" and the "dressed" fields. 

On the other hand, something along the lines of the second of the modifications, the in- 
troduction of the g- functional to achieve convergent expressions for the Feynman graphs, has 
been implicitly utilized previously. The reason for doing the calculation explicitly in terms of 
a g-functional rather than to regularize the divergences is to maintain control over the treat- 
ment so as to keep intact Noether's theorem and the other quantities and interrelations of the 
Lagrange-Hamilton based quantum theory. 

It will be demonstrated that in the present formulation the Lagrangian g-functional can be 
chosen such that without any regularization: 

(a) the Fock space is well-defined; 

(b) the propagators are of the Pauli-Villars form and are the Green's functions for the La- 
grangian g-functional constructed on the Fock space ; 

(c) the solution state vectors are given by a normalizable expansion on the Fock space basis; 

(d) the renormalization constants, as all the other quantities of the theory, are given unam- 
biguously by absolutely convergent Feynman - Schwinger integrals; 

(e) the solution fields, i.e., the fields which in the Heisenberg picture fulfill the equations of 
motion including the interactions, are non-local; 

(f) Haag's theorem @, §, || does not apply [§; 

(g) in the limit A — > (see below, Eq. ( |l.l| ) and the discussion) the renormalized S-matrix 
coincides with that computed in the conventional way; 

(h) the theory obeys all Ward-Takahashi identities derivable from Noether's theorem [|5|, ^, [j], 

1,1; 

(i) the concept of a wave function remains valid for relativistic QFT, and allows the description 
of many-body bound states; 

(j) the solutions are analytic functions of the coupling constant, say g 2 , as it tends towards 
zero; 
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(k) except for locality the theory fulfills the axioms of constructive field theory [10, 11], 



including that the physical quantities arising in QFT are at most tempered g-functions; and 



(1) Kailen's theorem 12] applies 



As is well known, the mathematical problems of the conventional formulation of QFT are the 
following: 

(A) The starting point are the solutions of the A = form (in contrast to the A — > limit) 
of the g- functional of our formulation. Most likely such solutions do not exist for non-trivial 
theories. Symptomatic for this are, for example, the statements: the S-matrix has an essential 
singularity at g 2 = 0; Fock space does not exist for non-trivial QFT. 

(B) In the regularization the Green's functions get mutilated by some prescription, be it the 
Pauli-Villars regularization or the dimensional regularization. Neither the regularized Green's 
functions, nor the remaining (finite) parts used in the calculation are the Green's functions of the 
problem at hand; they are disconnected from the assumed Lagrangian. 

(C) The Green's functions of the theory break the boundaries of tempered g-functions; as has 
been demonstrated earlier jj], ^] they contain meta-g- functions (see Appendix A). 

(D) According to the postulated properties of the (renormalized) fields, in particular their 
local character, Haag's theorem does apply and renders the results of QFT meaningless. 

(E) The problem of the convergence of the perturbation expansion of the S-matrix for a generic 
(3+1) QFT, the empirically observed at least semi-convergence of QED notwithstanding, remains 
unsolved and is not addressed in the present paper. 

In view in particular of the points (c) and (k) the mathematics of the present formulation is 
very much more simple than that of the conventional formulation. 

We now give a short overview of the points which will be demonstrated in the paper. 

The generalized (operator) functional, the g-functional, is defined by a sequence of functionals 
(see Appendix A): 

L{ipi, ...;A n } = L {ipi, . . . ; A„} + Li{i/) i} . . . ; X n } , (1.1) 

where A ra is a parameter (e.g., having the dimension of a length). The functional is constructed 
such that for n — > oo, i.e., for A n — > 0, L{X n } approaches the desired Lagrangian action (say, the 
quantum electrodynamics (QED) Lagrangian action), and that for A ^ all Feynman integrals 
of the Schwinger — Tomonaga expansion of the S-matrix converge. The "basic fields" t/ji etc., 
are quantized fields, i.e., operator fields, defined to obey given local commutation relations. They 
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can be chosen to fulfill the equations of motion resulting from Lq{-}. The (trivial) question of 
normalizability of the continuum can be taken care of in any way one likes, e.g., by the Weyl 
eigendifferential method, or by periodic boundary conditions. The basic fields can be used to 
construct a Fock space. No difficulty is encountered in the expansion of any quantities of the 
full theory in this Fock basis since, in view of the convergence of all graphs, all such quantities 
are mathematically fully defined directly in the Fock space without any mutilation by cut-offs 
or regularization of any integrals. The solutions, e.g., the expansion of the state vectors, can be 
obtained either by perturbative methods, or non-perturbatively, by matrix diagonalization. This 
way no difficulties arise in the description of bound states of composite systems, like, for example, 
positronium. 

In the present paper we provide a complete chain of the steps needed for the demonstration 
(of the mathematical aspects) of our formulation. However, we shall not give any proofs in terms 
of the rigorous mathematical language; our proofs, if at all given, will be in the simple language of 
a physicist. In fact, since the mathematical concepts arising in our formulation are rather simple, 
i.e., are at most tempered distributions (tempered g-functions), the need for the mathematician's 
language anyway is not that great. 

For concreteness we conduct the discussion in terms of the simplest possible non-trivial cases, 
i.e., a g- functional which has as its limiting form a spinor field interacting with a pseudoscalar 
field (Section 2), and the case of spinor electrodynamics (Section 3). The respective g-functionals 
are written in terms of Heisenberg-picture quantized fields. They are manifestly relativistically 
covariant in general, i.e., not only in the limit A n — > 0. To continue the development we need to 
derive the Hamiltonian, i.e., we must apply the variational calculus which in our case we need 
for quantized rather than for classical fields. As shown in Ref. JDI] the functional derivative, and 
hence the variational calculus, can be rigorously defined also for quantized fields, i.e., for operator 
fields and that then it yields the same expressions as for c-number fields. In view of this result, 
the conventional procedure which consists in splitting the Lagrangian into iff and , the 
"free" and the "interactions" part; in deriving the corresponding Hamiltonian; and in going over 
to the interaction picture, is mathematically fully defined and can be applied. At this point one 
can, as done conventionally, compute the S-matrix. All these steps are performed using only the 
basic fields. The solution fields, which are supposed to fulfill the complete equations of motion, 
are not required for computing the S-matrix. 

At this point one may terminate the development, as all perturbative predictions of the theory 
can now be computed. However, the physics and the need in particular for our modification (i) 
has remained hidden. This need becomes evident in the description of the solutions, in particular 
of bound many-body systems, as, for example, positronium. Thus, to demonstrate that the 
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modification (i) actually is basic to the viability of QFT, we continue the development, and we 
derive in Section 4 the form of the field operators, say ^, and of the state vectors, say \W), 
of a solution. In the g-functional formulation they are mathematically well-defined quantities. 
In particular, the solution fields can be expanded in terms of products of the basic fields; this 
demonstrates the existence of the Fock space for the interacting theory. Furthermore, this shows 
that the solution fields and the basic fields here are not unitarily equivalent; because of this Haag's 
theorem does not apply. The state vectors for any arbitrary physical system are normalizable; so 
are its wave functions, defined as usually as the matrix element (0| \W); this demonstrates the 
existence of the wave function in relativistic QFT. As in the case of the Pauli-Villars regularization 
procedure, for low energy processes the auxiliary (ghost, i.e., negative metric and/or "wrong" 
statistics, and normal) particle degrees of freedom for A n — > freeze out, except for providing 
for the convergence of all Feynman integrals. Owing to this convergence Kallen's theorem fl~2]| 
applies, which means that the A — > results are independent of the specific choice of the 
g-functional. 

As for previous work, as already mentioned above, the ingredient (ii) of the present approach, 
i.e., the use of modified Lagrangians, both polynomial and non-polynomial, to achieve convergent 
Feynman integrals is very old (l4|, [15], [l6|, . It has been employed in many papers; for example, in 
more recent work devoted to the proof of the existence of the cj) 4 theory [jl8| , |l9| 1 the development is 
carried out using regularized propagators without actually specifying the underlying Lagrangian; 
only the existence of such a Lagrangian is postulated. 

A discussion of the formulation of QFT in terms of path integrals lies outside of the frame of 
this paper. 

2 YUKAWA QUANTUM FIELD THEORY 

The simplest possible non-trivial quantum field theory is provided by interacting massive spin 
and spin ^ fields. This system has no further symmetries beyond those imposed by Lorentz 
invariance. The presence of further invariances, e.g., gauge invariance, requires imposing con- 
straints on the solutions. We shall treat such a case in the next Section. Here we consider the 
minimal theory, which, however, involves the full mathematical apparatus needed for the present 
re- formulation of QFT. 

Consider the Lagrangian action g-functional 

L{ij>...} = L F (X) + L B (X) + Lj(A) , (2.1) 
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(for simplicity we write A instead of A n ; also the limit A — > is understood throughout) 
Lp(X) = — d A x : {£i?/>i(x) (jd + mi) ipi(x) + £2^2(2;) (t<9 + m 2 ) ^2(2:) 



+ {jd + m 3 ) 4> 3 (x)} : (2.2) 

LbW = -J d A x : {i (d^x)) 2 + \ {M V {x)) 2 }: (2.3) 

U = : : . (2.4) 

This form is manifestly relativistically invariant. It contains 2 "auxiliary" Fermion fields. As we 
will see, this is a possible choice; there exist infinitely many other possible choices. 

The fields ipi(x), <p(x) are defined to obey the commutation relations = |e| = 1, see 
below) 

_L _L "1 q 

V> ia (x,t), typ{y,t) + e typi-y^t) V>ia(x,t)J = ZiSijSapS (x-y) (2.5) 

[^(x,t), ^(y,t)]_ = i<5 3 (x-y) , (2.6) 

where a, (3 are spinor indices. Fields which have £ = — 1 are ghost fields while normal fields have 
£ = +1, and e = — 1 yields the "wrong" statistics. Further, we define 

mi = m (2.7a) 

m 2 = y (2.7b) 

m 3 = y (2.7c) 

with C2 and C3 real positive constants. This defines the A — > limit of the g-functional. The 
c-number constants m will be determined below. 



Since, as shown in [|13|| , variational calculus is valid not only for c-number but also for operator 
fields one can go through the usual procedure, i.e., use the non-interacting part of L to obtain 
Hq, go over to the interaction picture, and compute the S-matrix in the standard manner (see 
Appendix B). To that end one needs the graph expansion of the Tomonaga - Schwinger equation. 
In view of the form of (|2.l| ), ([2.2|), this expansion will turn out to be exactly of the Pauli-Villars 
form as we now demonstrate. 

Consider the Fermion Green's function as it arises, for example, in the Wick expansion of the 
second order term of the Neumann series of the Tomonaga - Schwinger equation: 

G F = T «i«fc (o| M x ) My) |0> • (2.8) 

jk 
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In view of the commutation relations ( |2.5[) we have in momentum space (for brevity we omit the 
infinitesimal imaginary parts in the denominators) 

X 7P + mi jp + m 2 -yp + m 3 

Lrp = J^l — 9 + K 2 — 9 + A 3 



p 2 — m\ 



p 2 — m\ 



p 2 — m 2 



<yp UKt + K 2 + K 3 ) p 4 - [iTi(ml + ml) + K 2 (m 2 + m 2 3 ) + *C 3 (m? + mf) 
K\m\m\ + K 2 m\m\ + K 3 m\m\ j 



+ 



+ {(ifimi + -fT 2 "i2 + K 3 m 3 ) p 4 



+ 



+ 



K x m x (ml + ml) + K 2 m 2 (m 2 + m|) + i^ms^ + mi) 



K\m\m\m\ + K 2 m 2 m\m\ + K 3 m 3 m\m^ j 



(P -mfjip -m 2 ){p -mi) 



Here 



(2.9) 



(2.10) 



In order to cancel the terms of the numerator containing the factors (jp) p 4 and p 4 there must 
hold 



E Ki = 
^2 miKi = . 

Hence (we put K\ = 1 so as to achieve the needed form for A — ► 0) 

C3 — Ami 



K 3 



c 3 - c 2 
c 2 — Ami 



c 3 - c 2 

With these relations ( [2.9| ) can be simplified to the expression 



(2.11) 
(2.12) 

(2.13a) 
(2.13b) 



G F 



TP 



+ 



Ki ml + K 2 m 2 2 + K 3 m\ 



K\m\ + K 2 ml + K 3 m 3 



P 2 + 



P 2 + 



Ki K2 ^ K3 

2 ' 1 ' 2 

m\ 772-2 m 3 

K x K 2 K 3 

— + — + — 

m\ m 2 m 3 



222 
m l m 2 m 3 



222 
m x m 2 m 3 



x (p 2 - m?)(p 2 - m\){p 2 - mf} 
Here K\ has been retained so as to exhibit the symmetries of the expression. 



(2.14) 
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This way we have asymptotically Gf — ► \p\~ 3 ■ Together with the asymptotic character \p\~ 2 
of the Boson propagator one sees that the elementary Feynman graphs for the Fermion self-energy, 
the Boson self-energy, and the vertex correction, have asymptotic behavior as |p| -5 , H~ 6 , and 
|p|~ 8 , respectively. The integrals / d 4 p(-) thus are absolutely convergent, and no auxiliary (ghost) 
particles are required for the Boson field. The choice (2.2), ( p.3| ) for the g-functional thus turns 
out to be satisfactory. If one has chosen C2 < C3 then tp2 is a ghost field and £2 = — 1> £1 = £3 = 1 
in (ph. 



bi(p) ?4 a) (p) e^—*') + 4(p) v \ a \p) e -*(P*-***>l (2.15) 



The fields ipi(x) have the usual momentum space expansion: 



d s p Irrii 

(2^72 



,(«) 



where uii = + y p 2 + m 2 , and the operators b, d obey the usual, £- and e- modified anti-commutation 
relations 



^(p), ^ PJT (q) 



ZiSijSap 5 3 (p - q) 



(2.16a) 



and where all other anti-commutators vanish. 



Sij S a/3 5 3 (p - q) , 



(2.16b) 



One now can compute the renormalization constants, i.e., the mass shifts for the Fermion and 
the Boson, the vacuum-polarization shielding of the coupling constant, and the normalization 
constants, which all are finite, and proceed to compute the S-matrix. In view of the convergence of 



the renormalization constants here Kailen's theorem [12| applies; thus according to that theorem 
the renormalized results are independent of the renormalization scheme, which means that all 
choices of the g-functionals having the same limiting form and yielding convergent Feynman 
graphs will give identical results. 

We conclude this section by indicating the manner in which contact singularities arise in the 
Feynman - Schwinger integrals. To that end, consider the Fermion self-energy which is in second 
order 







ijkl 



x-y) S M (x - y)j5^j(y) 



(2.17) 



Here indicates that the integration has to be over an open domain so as not to include the 
contact terms H, |j, and ipi{x) and ipj(y) serve the role of test functions. The Wightman functions 
which retain the contact singularities are computed over the full domain and hence then is to 
be omitted. 
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To evaluate this expression with the aim of localizing the most singular terms we recall 



A(x) = JL {— log(-i/,V7 2 + *-s)+...j (2.18a) 
log(-z + ie) = log\z\ + iir9(z) (2.18b) 



,2 



Sij{x) — 5ijSj{x)Kj 



(2.19a) 



Sj(x) 



2jx 



+ 



47T 2 1 (^ 2 — *e) 2 2 x 2 — ie 



7X 



+ 



\ log ( — i to 2 x 2 7 2 + ie 



, , , , ( 2 - 19b ) 

X^ — IE 

Here 7 = 1.781 ... . Inserting the expansions (|2~T^ ), (gl9b|) in (|2~17|) we obtain, putting y = 
in view of translational invariance 



| d 4 xi;(x) l5 J2 K J 



2jx 



+ 



1 



— m 



7X 



(x 2 — ie) 3 2 (x 2 — ie) 2 
"fx ( 1 



J (x 2 — ie) 2 4 (x 2 — ie) 2 
2 



log 



2 2 ~2 1 • 

/i x 7 + ze 



1 



~r 2 " ■ lo § 

4 x z — ze 



j m 2 x 2 7 2 + ie ] + . . . [>75 ^(0) 



(2.20) 



Prom (gig) and ( gig) we see that with the choice of the constants Kj of ( [2.10D the first two 
terms of the "individual" propagators Sj(x) cancel in the summation of ( [2.20D and do not appear 
in E. However, in view of our present task we shall look at one of the propagators, Sj(x), leaving 
the cancellation to the end. 



Following Ref. [|J we define x^{ x ) such that there holds 

4>{x) = x»(x) + -0(0) 



(2.21a) 



x^{x)+ I d A x'^(x')6\x') 



(2.21b) 



Consider now the second term of the left-hand side of ( p. 20 ). Inserting (2.21a) in that term we 
have 



T. ~ 

1 3 — 



d 4 x Xfi(x) 



(x 2 — irj) 2 
d 4 x' 4>(x')5 4 (x')i;(0) 



(x 2 — irj)* 



(2.22) 
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where we have used the version fl2.21b| ) in order to make manifest the contact character of this 
term. The second term of ( [2.22 ) shows that /, and hence the individual terms of S, can contain 
contact singularities with amplitude proportional to 

1 



d x 



(2.23) 



(x 2 — irj) 2 

As shown in detail in Ref. [0, Eq. ( |2.23| ) indeed is logarithmically divergent. However, being a 
contact term it can contribute only to the Wightman function, and not to the S-matrix. At any 
rate, here it is eliminated upon summation over j in ( |2.20| ), also for the Wightman function. All 
remaining terms are regular for x^ = and are best evaluated in momentum space. Their sum 
yields convergent integrals, individually for each j. 

To summarize this Section, we have seen that in the g-functional formulation all matrix 
elements of the S-matrix are given by convergent expressions and hence are mathematically 
rigorously defined, and the Green's functions are those associated with the Lagrangian. In this the 
present formulation differs from the conventional regularization procedure where the regularized 
Green's function is disconnected from the Lagrangian and where the choice of the regularizing 
function is arbitrary and only dictated by convenience. At any rate, in view of Kallen's theorem 
[12 1 the results of the conventional procedure turn out to be "correct" in the sense that they agree 
with those obtained as rigorous solutions of the g-functional formulation. 



3 QUANTUM ELECTRODYNAMICS 



As another example we now consider QED. Thus we have to take care of charge conservation, in 
addition to Lorentz invariance. We take here the g-functional 



with 



L F (X) 
L B 
Li 



and 



L{V---} = L F (X) + L B (A) + L/(A) 



d 4 x : [ipi(x)(jd + mi)ip(x) - ip2(x){jd + m 2 )V>2(£)] : 



J d 4 x : X - F^(x)F^{x) - £ (d, L B v d^B u - M 2 B v B l 



i4 / d 4 x : [(ei^n^i) + (e 2 ip 2 lM] {A^+B^) : 



(3.1) 

(3.2) 
(3.3) 
(3.4) 

(3.5) 
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This g-functional thus contains besides the Fermion and the photon field two ghost fields: one 
Fermion and one Boson. The latter is a mixture of spin and spin 1. 

We check for charge conservation. To that end we introduce the substitution 

ih -> e ia fa , ^ 2 -> e lf3 V>2 , (3.6) 

with a, (3 infinitesimal, and we demand 

SL 5L . . 

e - ■ Tp = • (3 - 7) 

From Noether's theorem we have 

d M = , d, = , (3.8) 

ffi = * e n $n Ifi *Pn , (3.9) 

i.e., separately conserved currents. With the ansatz Eq. Q3.4| ) here, in contrast to the Pauli-Villars 
procedure, no photo-excitation from ipx to ip2 exists. Therefore one is free to chose the ghost 
coupling constant, e2 at will. Thus in the Fermion mass-renormalization graph asymptotically at 
large loop momenta the Fermion propagator retains the character | p | . 

The only place where the Fermion ghost contributes is in the vacuum polarization, in that 
for every closed Fermion particle loop there exists an analogous ghost loop. The lowest order 
contribution then is 



n = 

where 



r 2 



d*p (3.10) 



q = p-k (3.11) 



and where the minus sign of the ghost contribution results from the ghost field having the "wrong" 
statistics, i.e., e = — 1 in ( |2.5| ). In each of these terms the quadratic divergence is a contact 
singularity and does not contribute to the open integral, while the surviving terms are 

e 2 / d 4 p PWi-C 2 m|) + ••• (3 12) 

J (p 2 — m\)(q 2 — mf)(p 2 — m?,)^ 2 — m?,) 

Choosing 

C 2 = 4 = 4 (3-13) 
ef 

the integral is absolutely convergent. With this choice in the limit m 2 — ► oo the Fermion ghost 
not only freezes out but also decouples. 
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We now turn to the Boson sector. The photon propagator is the usual one (Feynman gauge) 

Ap h = / d 'P ( 3 - 14 ) 

J p — irj 

while for the ghost it is directly the Klein-Gordon propagator 

S. 



A 



G 



The full (photon + ghost) Boson propagator thus is 

j\ — _ f *W ^ 2 (3 IQ'j 

J p 2 (p 2 — M 2 ) — ir\ ' 

this way together with the Fermion propagator (j) — m) , all graphs are absolutely convergent 
for \p\ — > oo, and the limits m 2 , M 2 — > oo can be safely carried out. 



4 FORM OF THE STATIONARY STATE SOLUTIONS 



Having seen that the S-matrix can be computed in a mathematically solid manner, we now 
shall show that the solutions themselves are well defined and can be used to construct the wave 
functions of the system, which then can be used to evaluate its characteristics exactly as in the 
familiar non-relativistic quantum mechanics. It is here that the modification (i) mentioned in the 
Introduction becomes explicitly evident. 

The computation of the S-matrix in the interaction picture is based on a time-dependent 
treatment. Of course, this treatment can be applied to the computation of stationary states 
. On the other hand a stationary state is easiest to visualize in the time-independent 
treatment. Therefore, we will begin by illustrating the form of the solution in the Schrodinger 
picture using as examples the case of a single (physical) particle of the Yukawa theory given 
above, and that of a composite bound state system for the case of different Fermions, for example 
the bound state of a muon and an electron. We then will re-state the problem in the interaction 
picture and derive a covariant eigenvalue equation in terms of the U-matrix elements. Of course, 
for a convergent (or, within its range of convergence, for a semi-convergent) theory the results 
are independent of the choice of the method. In particular, it is well known that the Schrodinger- 
picture treatment is fully equivalent to the treatment by the Feynman graph expansion, only 
being enormously more cumbersome for performing an actual calculation. 

In the Schrodinger picture the state vector of a given system, |VF( s )(t)), obeys the equation 

1 \wW(t)) = -iH(\) \W( s \t)) (4.1) 
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where 

H(X) = H (X) + H 1 (X) 
is the Hamiltonian of the system; Hq(X) and -Hi (A) are associated with Lq(X) 



(4.2) 

L F (X) + L B (X) 

and L/(A), of Eqs ( p.l| ) through (2.4), respectively. The basic fields are those of the Heisenberg 

picture taken at the time t = 0. They are taken to be solutions of Hq(X). So, for example, the 

(s) 

Fourier decomposition of the Schrodinger-picture fields ip\ (x) is 



4 s \x) 



4 +){s) w + 4~ ){s) (x) 



d 3 p 




br,u D e ipx + dtv D e- ipx 



(4.3) 



(2vr)3/2 V Ui V p " p ~ 1 "p^ p 

where we have suppressed the spin indices, and similarly for the other basic fields of the La- 
grangian. Also, we use the notation uj{ = + \Jp 2 + fnf- In ( |4.3| ) we have introduced the usual 
notation for the positive frequency {^ + \ annihilation operator) and negative frequency (ip^~\ 
creation operator) fields. In principle the basis ( fl.3| ) must be discretized in some manner, e.g., 
by Weyl's eigendifferential method. From now on we drop the superscript (S) indicating the 
Schrodinger picture, and the explicit reference to the parameter A. 



We note the commutation relations 
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Z5 a(3 <5 3 (p-q) 



(4.4) 



dap, rfjjq 



ZSaf, <5 3 (p-q) 



(4.5) 



where as above £ 
statistics. 



= £<W<5 3 (p-q) 
+ 1 for particles and = — 1 for ghosts, and e 



(4.6) 

1 for the "wrong" 



A solution for a stationary state of a given system can be obtained by expansion in the 
Fock-space representation, i.e., by constructing a complete basis of states having the quantum 
numbers of the desired state (momentum, spin, etc.) and diagonalizing the total Hamiltonian, 
Eq. (f4.2j), in that basis. The Schrodinger-picture state vector for the system consisting of a single 
Fermion of momentum p at infinite distance of any other particles then will have the general 
form (suppressing the spin indices) 



\W(t)) 



-iEt 



\W) 



(4.7) 



\W) 



C wo b P |0) + J d 3 k doi(k) bl_ k 4 |0> 
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+ J d 3 h d 3 k 2 C 102 (kik 2 ) bl_ ki 4 x _ k2 4 2 |0> + . . . 
+ J d 3 Pl d 3 P2 C 2 i (piP 2 ) & P _ Pl 6^ 4,2 |0) + . . . 



+ ... (4.8) 

where actually, e.g., 6 k = £^ &j k to account for the auxiliary fields as implied by (|2.2|) . In (|4.7|) , 
-E 1 is the energy, i.e., E = +\/p 2 + M 2 , where M is the at this juncture unknown mass of the 
system. 

For the composite system we assume two distinct Fermions but only one Boson species, and 
their assorted auxiliary particles. The state vector for the composite system is similarly to above 



\W&) = j d 3 PiCnooo(pi) b^ pi bW |0> 

+ J d 3 Pl d 3 k Cnooi (pik) b^ pi b p % 4 |0) + . . . 

+ J d 3 Pl d 3 P2 d 3 P3 C 2211 o(piP 2 p 3 ) ^4 b p b ll p2 6gi p , 4f |0) + . . . 

+ ... (4.9) 



In the usual manner we now compute the Hamiltonian matrix by inserting (4.8) in (4.1) and 
by multiplying on the left with the Hermitian conjugate of one of the components of the state 
vector at a time, for example, {{f3}\ = (0| 6 p _k a k- This way we obtain the matrix equation 

(E-E m ) C m = E / d'h ... C {a} ({/3}\ Hi \{a}) (4.10) 

{a} 

where the fact has been used that Hq is diagonal in the Fock space configurations, yielding the 
value E{py for the configuration {/?}. 



The matrix (p~Tc| ) now can be diagonalized, either directly, i.e., non-perturbatively, or pertur- 



batively, say, by the Rayleigh - Schrodinger method, as a power series in the coupling constant, 
say g, i.e., by writing 

E = E + gE 1 +g 2 E 2 +g 3 E 3 + ... , (4.11) 

where for the case of semi-convergence the series has to be terminated. (Equivalently, in the 
non-perturbative treatment the Hamiltonian matrix has to be truncated.) One has, for example, 

2p V^r x [,r 1 ({l}|ffl|{Q» ({a}jffl[{l}) 

g E 2 = 2_^{a\ J d{a} _^ (4.12) 
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where {1} is the ground configuration, i.e., the first term of (4.8). The integration in ( [4.12 ) 



is over the quantum numbers of the configuration {a}, for example, the momentum k for the 

b ] a ] - 
5 P -k "k- 



configuration {a} = &I,_ k aj< ' 



d{a} ({b j P }\ Hi |W) ({«}| Hi \{bl}) = r 3k ({^}|^i|^_ k 4)(^-k a kl^il{&p}) 



E a — E J LUk + Wp-k — E 

(4.13) 

In the conventional formulation integrals of this type diverge. Owing to the auxiliary parti- 
cles they are non-divergent in the g-functional formulation; thus Fock space is well-defined. As 
long as the configurations {a}, {/?}, retained in the Hamiltonian matrix belong to those appear- 
ing as intermediate states in the Feynman graphs connected with the terms of the convergent 
part of the presumably semi-convergent series expansion for the S-matrix, they have physical 
significance. If they belong to that part of the series expansion where the terms have started to 
grow, their retention would tend to diminish the accuracy of the results; and thus they should 
be omitted. This condition thus represents the criterion for the truncation of the secular equa- 
tion, Eq. (4.10), needed for the direct diagonalization. (Further details of the solution of field 



theory problems by diagonalization in the Schrodinger picture, in particular the treatment of the 



relativistic center-of-mass, have been presented in Refs. 24, E5[. 



Since the Feynman integrals converge and themselves are independent of g 2 , E2 is independent 
of g 2 , and E n is independent of g n . Hence, the eigenvalue E of the solution joins smoothly to the 
non-interacting theory. More precisely: if in fact the S-matrix expansion is only semi-convergent, 
then as g 2 becomes smaller the beginning of the growth of the terms in the expansion gets 
postponed, and for g 2 — > this beginning occurs at arbitrarily large order n. Thus, as mentioned 
in point (j) in the Introduction, E, Eq. ( 4.11| ), and also the S-matrix itself, in the g-functional 



formulation is analytic at the point g = 0. 

These analytic properties of the point g 2 = are seen by the following limiting procedure: 
Given a value of A, and a value of N, then in 

E N ' X (g) = E 9 n E x n (4.14) 

there exists a go ^ such that for \g\ < \go\ and n < N there holds 

n\g n E x n < (n-k)lg n - k E x n _ k (4.15) 

for positive integer k. Thus for g — > the value of E N '^(g) and its n-th derivatives with n < N 
are independent of N. 

We note the form of the wave function of the system. It is given by [re-introducing the 
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superscript (£)] 

v (xi,t) = (0| V {s) \W {s \t)) 
where for the single Fermion case ty^) [ s field 



(4.16) 



IK 



JK 



*(*> = e n ^ {+) (^) n ^ (+) (x,) n 

{k} i j k 



(4.17) 



For the composite system the field is 

ift Jre «k Z« "Ik 

*w = y: n ^ {a)(+) (^) n ^ (a)(+) (x;) n ^ )(+) (**) n n ^ (+) ( x -) • 

{k} i j k I m 

(4.18) 

For the one-Fermion system only the following terms of the full set of configurations, are 
needed 



(S) 

(1) 



V (+) ( Xl ) [l + ^(+)( yi ) + ^ +) (yiV (+) (y 2 ) + . . . 
+ ^(+)( X1 ) v (+) (x 2 ) V (+) (x 3 ) [i + ^ (+) (yi) + v (+) (yi) ^ (+) (y 2 ) + ... j 

+ ... . (4.19) 
Analogously only a subset of the configurations of Eq. ( |4.18| ) is needed for the composite system. 



According to ( j4.7| ) and ( |4.8[) the wave function rj(xi,t) of ( f4.16| ) has the overall time dependence 
e -iEt anc j j g a m ixt ure f configurations containing different numbers of quanta of the basic fields. 
This is fully analogous to the configuration mixing in atomic or nuclear physics. One may call 
ty( s } of (4-17) the general Schrodinger-picture solution field, and \W^ s \t)}, ( |4.8j ), the Schrodinger- 
picture solution state vector for the one-Fermion system. And, given that the basic fields obey 



local commutation relations it is obvious that this is not the case for the solution fields (4.17) or 



- which is the equal- 



( 4.19| ). As can easily be verified, the commutation relation ty( s \ 
time relation of the corresponding interaction picture fields - is a mixture of c-numbers, creation 
operators, annihilation operators, and number operators, and is non-local; i.e., these fields do not 
commute at space-like separation, which is contrary to one of the axioms of constructive field 
theory. This way one of the basic assumptions of Haag's theorem, [jl], ||, ^] i.e., that the equal-time 
commutation relations of the fields Vl/ yield a local c-number, is not fulfilled; hence, that theorem 
does not apply in our g- functional formulation of QFT, and hence the interaction picture exists 
for non-trivial theories. 

As the last point we now indicate that the wavefunction can be computed also in a covariant 



manner. To achieve this we now have to re-express (4.1) in terms of the U-matrix. This will 
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also show the connections with the results of the previous Sections. To that end we re-write 
our expressions in the interaction picture. We define for the state vector of the system in the 
stationary state n of energy E n 



|W«(*)> = e- iE ^ \W<P) = E ir 

{a} 



WOO, 



(4.20) 



Here we have compactified the notation of (fOj) by introducing the combined index {a}. The sum 



over {a} in (4.20) thus contains both the summation over the configurations and the integration 
(actually summation for a Weyl representation) over the momenta. The product of the creation 
operators of ( [4.8| ) is denoted by Similarly, we re- write eq ( 4.17| ) for the interaction-picture 

fields as 



{a} 



= E II - E*£1 (I W) 



{a} abc 



{a} 



(4.21) 



taken at the same time t for all fields. In (4.21) the indices, a,b,c, denote also the momenta in 
the field expansions according to (4.3), for example. 



Herewith, we have 

TfcOM) = <0| |W«) 
Being in the interaction picture the time-dependence of the products of fields is 



(+)(!) 
M 



(4.22) 



(4.23a) 



E< 



{a} 



■ J{a} = ^ *i (4.23b) 

i 

where the Ei are the energies of the individual fields participating in the configuration {a}, 



eq (4.21). Hence, the amplitudes of the state vector for our state must have the time-dependence 



,WC0, 

J {a} 



-i{E n ~E {a} )t («)(/) 
{a} 



(4.24) 



in order to achieve the overall time dependence of ( 4.20j ). 



We now follow the development of the state beginning at time to by means of the evolution 
operator U. Dropping the superscript (I) we have 



|W n (ti)> = U(t u t ) \W n (to)) 



or, in full detail 



P 



(4.25) 



(4.26) 



17 



Combining ( [4.26 ) with ( 4.24j) we obtain 



e -i(£„-E {a} )(ti-to) w (n) 



,(n) 



(4.27) 



dropping the argument to m the amplitudes. We now introduce the operator 



e 'o 



-i f' 1 (£„-#<,)<**' 



U<®(E n ;ti,to) 



(4.28) 



which is diagonal in our basis. It has the structure of a mass counter term. Herewith, Q4.26 ) can 



be written as 



E [ s °fi U ap( E n, h - t ) - U aP {h - t )] w$ } = (4.29) 



n 



which is a (non-linear) eigenvalue equation for the energy and the state vector amplitudes. Upon 
the limit to ~* —oo, t\ — > +oo, one arrives at the covariant expressions for the S-matrix elements 
which can be computed by the familiar Feynman graph expansion. This way we see that the 
concept of the wavefunction of a system is not limited to non-relativistic physics but has a precise 
meaning also in the frame of relativistic g- functional QFT. 



In this paper we have demonstrated that in the g-functional formulation of quantum field theory 
every step leading from the Lagrangian to the S-matrix can be defined unambiguously in a 
mathematically rigorous way; the renormalizations are non-divergent and can be carried through 
explicitly; and the state vectors are normalizable. As mentioned in point (E) of the Introduction, 
concerning the convergence properties of QFT we assume semi-convergence of the perturbation 
expansion. The solution fields are non-local and hence Haag's theorem does not apply; i.e., the 
interaction picture exists. The state vectors of the solutions are linear superpositions of Fock space 
configurations. This is true both for single particles, e.g., electrons, and for composite systems, 
e.g., positronium. For a small coupling constant, in the solution the ground configuration has the 
largest amplitude; the higher configurations are corrections. The solutions thus connect smoothly 
with the "free" theory; i.e., for a fixed A the limit g 1 — > within the limitations imposed by the 
possible semi-convergence of the expansion is analytic in the usual sense: for given e and n, 
there is a g 2 such that for g 2 < g 2 , the n th configuration has an admixture \A n \ 2 < e. Herewith 
we have shown that non-trivial relativistic QFT is not only consistent mathematically, but also 
makes sense from the point of view of physics. Thus, the theory is "transparent," i.e., from 
contemplating the solutions one can draw conclusions concerning the underlying Lagrangian. 
This connection is not colored by the mathematics: the result does not depend on the choice 
of the g-functional, as long as it has the desired Lagrangian as the limiting form, and it has 
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the needed convergence properties. Conversely, one then is on safe grounds to judge whether a 
chosen Lagrangian, i.e., always defined as the limit of a g-functional, is suitable for describing the 
physics at hand. Thus we can conclude that non-trivial relativistic quantum field theory exists 
as a rigorous mathematical theory and that it can be employed in attempting the description of 
Nature. 



APPENDIX A 



g-functions and meta-g-functions 

Generalized functions (g-functions), also called distributions, are defined in a limit procedure 
by their action on a function, the "test function." In the terminology of Lighthill |2(| the limit 
procedure is carried out in terms of sets of "good functions" which contain a parameter, A, and 
the g-functions arise in the limit A — > 0. Considering as an example the 5— function, we have 

Six) = MmAxix) (A.l) 

A^O 

and 

Jim f A x (x)f(x)dx = /(0) . (A.2) 

A — >0 J 

Of course, ( |A.l| ) is only a symbolic notation and to be understood in the context of ( |A.2| ). There 
exist unlimitedly many different choices for the sets of "good functions" A^(x) which can be used. 
But all of them make mathematical sense only for A / 0. The functions A\(x) are meaningless 
for A = 0. The limit procedure thus is to be carried over the domain [Ao, 0) for A, i.e., the domain 
open at A = 0, and where Ao is some arbitrary, "small" value. Also, given a set A\(x), the 
functions f(x), the "test functions," must have certain characteristics for (|A.2|) to hold. 



All this is well-known; the mathematics of g-functions is essentially complete p6|, p7]]. The 



most essential of their characteristics is that their action, analogous to (A.2), is unambiguous, 
and that they have unambiguous Fourier transformations. Thus the Fourier transform of (A.2) 
is given by the convolution of the Fourier transforms of A\(x) and of f(x). 

When attempting to apply (|A.2j) to the case 



x \x 
[g(x) a "good" function] one arrives formally at 



f(x) = - g(x) ~ (- + Z5(x) ) g{x) (A.3) 



lim J A\(x) — g(x) dx = lim J A\(x) ( — + Z 5(x)^j g(x)dx (A. 4) 
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which breaks the frame of g-functions; it involves the product of two g-functions. Such quantities 
have been called "meta-g- functions." [j/J The Wightman functions of QFT frequently contain 
meta- g-functions. 



Going back to (A.l) one finds 

[5(x)] 2 = hm A Al (x)A A2 (x) . (A.5) 



A 2 ^0 



One has here two independent limits. The result one would obtain from 

lim f A Xl (x) A A2 (x) g{x) dx (A.6) 

Ai,A2— >0 J 



depends on the manner one performs the limit procedure. Thus (A.6) is inherently ambiguous. 
Meta-g-functions unavoidably yield ambiguous results. Not much is known about their prop- 
erties. Still, one can define their "singularity character," which is independent of the way one 
performs the limits, e.g., in ( |A.6| ). Thus one may call a meta-g-function D^ n \x) to be of degree 
n (n =integer) if there holds 

J D^ n \x) x n+k dx = (A.7) 
with integer k > 0. With ( A.7 ) one sees that the general solution of the equation (measure dx) 

x n f(x) = g(x) (A.8) 

in addition to g-functions (5— functions and their derivatives) contains also meta-g-functions of 
degree up to n. (The Eq. ( |A.8| ) can be used to define meta-g-functions of degree n Q]. An example 
of a degree 2 meta-g-function is Eq. ( |A.6| ). Meta-g-functions of degree n = 1 are g-functions.) 
Some further discussion on the properties of meta-g-functions is contained in Ref. M. 

Contact type g-functions and meta-g-functions. 

The open integral, f% dx-, is defined to leave out the single point x = 0. Precisely, it is 
defined as a definite integral over the range 0(a, b) = (a, 0),(0, b). It differs from the integral 
J dx ■ which has the range (a, b) = (a, 0] (0,6) = (a,0) [0,6). In order for these two integrals to 
be different, the integrand must be a (meta-)g- function of the contact type, for example, S(x), or 
[5(x)] 2 , etc. Thus, for example, 

f dx [/(x) 6{x) + g(x)] = f dx g(x) (A.9) 



for f(x) and g(x) (ordinary) test-functions. The g-function 8'(x) is not of the contact type as 
can be seen from the definition 

6>(x) = lim 6(X + V)-S(x-Tl) {A W) 
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It is retained in the open integral. Further details are given in Ref. J?j. 

We could carry out our analysis of the meta-g-functions because we used the procedure of 



Lighthill [26|. The Schwartz - Gel'fand definitions [27] do not lend themselves to the required 
generalizations. 



g-functionals and meta-g-functionals 

In analogy to the g-functions one can define generalized functionals, for short g-functionals. 
They are functionals of the fields ip(x), which can be c-numbers or operators. In analogy to the 
good functions generating the g-functions, the functionals generating the g-functionals contain 
a parameter A, and they are well-defined for A / 0; that means that the functional equations 
they are associated with (e.g., the equations of motion) have solutions for tp (and for the state 
vectors upon which the field operators act) for A / 0. They make no mathematical sense for 
A = 0. The g-functionals are then defined in analogy with the g-functions: in order to achieve 
mathematically sound results the calculation is performed for a non-zero, fixed value of A. Only 
after the completion of the calculation one goes to the limit A — > 0, in the strict sense that A = 
is not permitted; the limit procedure is to be carried out over the domain [Ao,0), which is open 
at A = 0, with arbitrary "small" Ao- To emphasize: the limit procedure is carried out for the 
solutions of the theory (e.g., the S-matrix), i.e., the solutions obtained for "small" A. 

Furthermore, in complete analogy to g-function theory one can give for typical, general cases 
of g-functionals rules, for example, Feynman rules, on how to write the results of the A — > limit 
procedure. And, in the same way as for g-functions, one always can check the rules by actually 
computing the A — ► results. 

As is well-known, the solutions of the g-functional theories contain g-functions. For example, 
the Feynman integrals, say for QED, contain light-cone 5— functions and their derivatives. These 
g-functions survive the A — > limit procedure. One may say that g-functionals lead to g-functions. 
Inasmuch as g-functions have well-defined characteristics, the solutions of g-functional theories 
are also fully, unambiguously defined. 

One also can introduce the concept "meta-generalized-functionals" arising in the context of 
g-functionals as analogues of the meta-g-functions which arise in the context of g-functions. In the 
same way that meta-g-functions are inherently ambiguous, the meta-g-functionals are inherently 
ambiguous. Examples for the meta-g-functionals are the "non-renormalizable" quantum field the- 
ories: they cannot be made unambiguous - even by a A — > limit procedure. In fact, the Feynman 
integrals arising in non-renormalizable theories contain non-contact-type meta-g-functions, which 
belong in the solutions Q and which render the solutions ambiguous. Again, one may say that 
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meta-g- functional theories lead to meta-g- functions. 



APPENDIX B 



Negative Metric Fields 

Consider the action 

L = j ip(x) (jd + m) ip(x) d 4 x 
The Fourier expansion of the fields is 



(B.l) 



ijj(x) 



m 



d 3 p 

WW 2 



b p u{p) e ^ px - Et ^ + dj, v(p) e -i{^-Et) 



(B.2) 



The Fock operators obey the usual anti-commutation relation except for the "wrong" sign: 



The Hamiltonian here is 

H 



The equations of motion are, 

d 3 p f d 3 k 



S i(px.—E p t) 



We have 



: b^bp : 6 k - fe k : fe p fe p : = : b p b p : 6 k - (-6 k & P - fe p 6 k + 6 p 6 k ) 6, 



: 6 p 6 p : 6 k - 5 3 {k - p)6 k + 6 p 6 k 6 P 



: 6 P 6 P : b k - S 3 {k - p)b k - : 6 p 6 p : b k 



= -,5 3 (k-p)6 k 



and similarly for the other terms. Thus 

d 3 p 



[H,1>(x)\- = -Jj^E p {b p u p e^- E ^-dlv p e-^- E ^} 



1 d_ 

~i di 



il>(x) 



(B.3) 



(B.4) 



PM*)}- = -Jj^Jj^ E A[- b i b »-- b *- b *-- b i> b »\ 

: d\,d p : 4 - 4 : d) p d p :] Wpe^^')} . (B.5) 



(B.6) 



(B.7) 



which agrees with the equations of motion for positive metric fields. The expressions for the 
interaction picture development thus remain unchanged. 
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